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Abstract 

The static potential between an infinitely heavy quark and anti- 
quark is derived in the framework of perturbative QCD to three loops 
by performing a full calculation of the two-loop diagrams and using 
the renormalization group. The contribution of massless fermions is 
included. 

The force law between infinitely heavy quarks has been investigated since 
more than 20 years because of its importance for a deeper understanding 
of the strong interactions. The static quark-antiquark potential is a very 
fundamental concept, constituting the non-abelian analog to the Coulomb 
potential of electrodynamics, and also enters as a vital ingredient in the de- 
scription of non-relativistic bound states like quarkonia. It is widely believed 
to consist of two parts: a Coulombic term at short distances which can be 
derived from field theory by using perturbative QCD, and a long-ranged con- 
fining term whose derivation from first principles presumably requires much 
more advanced methods. Although an analysis based on perturbation theory 
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alone thus cannot give the complete potential, the result of such an effort 
would nevertheless be very useful. It could provide an improved input for 
QCD inspired potential models or even describe very heavy systems to a rea- 
sonable accuracy by itself. It could be compared with the potential obtained 
from numerical studies using lattice gauge theory and it might also give some 
hints on the nonperturbative regime. 

The first investigation of the static perturbative QCD potential has been 
performed in Although this work has been extended by several groups 
shortly after [§, [|, [3], [5], and some of these also studied aspects of the two- 
[0] and even three-loop diagrams Q, there is still no full calculation of the 
two-loop diagrams available. The purpose of the present paper is to fill this 
gap and hence, by exploiting the renormalization group equation, to obtain 
the three-loop potential. 

Before turning to the actual analysis, let us first recall the calculational 
procedure employed. It seems appropriate to begin with the simplest case, 
the Abelian theory without massless fermions. 

The static potential in QED can be defined in a way which makes its gauge 
invariance manifest via the vacuum expectation value of a Wilson loop taken 
about a rectangle of width R and length T ^> R: 



where P denotes the path ordering prescription. 

The functional integral can be calculated exactly, and one indeed finds the 
Coulomb potential plus an additional term which represents the self-energy 
of the sources 0. To compare with the non- Abelian theory it is, however, 
useful to go through the perturbative analysis as well. The Feynman rules 
for the source are as follows: a source-photon vertex corresponds to a factor 
iev^, with an additional minus sign for the "anti-source", and the source 
"propagator" reads 



in coordinate space or 



in momentum space. The four vector v is given by = (1, 0) and has only 
been introduced for notational reasons. The appearance of a propagator for 
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Figure 1: One-loop ladder diagrams. The double lines indicate static sources. 

the sources is a consequence of the time ordering prescription in the path 
integral, which introduces G-functions when expanding the exponential. 

Looking at the problem in this way shows the connection to another ap- 
proach: the potential beyond the infinite mass limit or for non-singlet sources 
is frequently also derived from the scattering operator; see for example |5], I7f]. 
The static QED potential thus should be derivable from the scattering op- 
erator of heavy electron effective theory, the QED analog of heavy quark 
effective theory, and the Feynman rules given above obviously support this 
point of view. 

Some care is required because the Feynman diagrams do not directly 
correspond to the potential but to exp[—iV(R)T]. The consequence is that 
in the abelian theory the one-photon exchange amplitude already gives the 
final result: 



-ie 2 Jdx dy v ll v v D tlv {x - y) — -ie 2 T J 



d 3 q e 



tqR 



(2vr) 3 -q 2 



where D^ u is used to represent the photon propagator and x = 0, y = R are 
understood. 

At one loop order, one encounters self-energy and vertex corrections which 
cancel due to the Ward identity, and the ladder diagrams shown in Fig. |T] 
which are best analyzed in coordinate space. Because of the simple struc- 
ture of the source propagator, Eq. (g), only integrations over time- variables 
remain. Adding the two diagrams removes the G-function corresponding 
to the anti-source propagator, and adding them once more with x *->■ x', the 
source propagator can also be removed and the one-photon exchange squared 
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is obtained, 

2 x Fig. | = ( - ie 2 T J dx dy D 00 (x - y)f . 

This behaviour of the ladder diagrams persists in higher orders 0, the ex- 
ponential thus starts to build up. 

To see the exponentiation in momentum space is more difficult, and re- 
quires that we specify the gauge — which will be Feynman gauge — and 
the special kinematic situation. As the sources are infinitely heavy they may 
carry any three-momenta without moving, but the actual values of these 
three-momenta are irrelevant as the only quantities that enter the calcula- 
tion are the momentum transfer q and the energies of the sources. The latter 
are required to vanish by the on-shell condition (implied in the propagator, 
for example) and consequently the only dimensionful parameter that remains 
is q 2 = — q 2 . Using dimensional regularization with D = 4 — 2e to handle 
infrared divergencies, the individual amplitudes for the diagrams in Fig. [1 
thus read 

Fig. 1(a) = ie*(±y±T(e)0yr(v.v'), (4) 
Fig. = ie ^±y± 2m (^yr(-vv% (5) 

with 

r(w) = In (w + V w 2 — l) , w > 1, (6) 

\Jw 2 — 1 v ' 

Z7T 

r(w) = —r(—w) H — , w < —1. (7) 

Vw — 1 

v 1 is used to denote the four-velocity of the anti-source, hence we should 
take the limit v' — > —v, v ■ v' — > —1. But this results in a badly diverging 
imaginary part for the uncrossed diagram. Keeping a relative motion by 
setting v' = (— yl + u 2 , u) we can, however, recognize this divergence as 
resulting from the Coulomb phase ||. Another way to see this fact would be 
to keep the kinetic energy in the heavy electron propagator. As the real parts 
of the diagrams cancel we thus again find that they are merely an iteration 
of the one-photon exchange. A similar analysis should be possible for the 
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higher order ladder diagrams of course, but it is obvious that the coordinate 
space approach is much easier in this respect. 

The inclusion of rif massless (i.e., m 2 <C q 2 ) fermions, although it makes 
an exact solution impossible, presents no problem in perturbation theory. 
The fermions appear as loops in the photon propagator and induce light-by- 
light scattering and in this way lead to an effective running coupling constant, 
i.e., 



Note that this effective coupling differs from the usual running coupling in 
the MS-scheme. Light-by-light scattering, in fact, first enters in three-loop 
graphs and is thus beyond the scope of this paper. 

When turning to the non-Abelian case, the Wilson loop must be gener- 
alized to 



where the matrices T a denote the group generators. Consequently the po- 
tential for a quark- ant iquark pair in a color-singlet state can be defined as 



In principle there are some problems connected to this definition, caused 
by the non-trivial topological structure of non-Abelian theories, which are, 
however, absent in the purely perturbative approach. 

As there is no way known to solve the QCD functional integral exactly, 
one has to resort to a perturbative treatment, which is, of course, more com- 
plicated than in the Abelian case: additional diagrams appear due to the 
trilinear and quartic gluon self couplings, and the presence of the generators 
in the source-gluon vertex influences the exponentiation as will be demon- 
strated. 

We will use Feynman gauge and the kinematics as described above again. 
Because the individual loop diagrams contain both infrared and ultraviolet 
divergencies, dimensional regularization will be employed, without, however, 
explicitly distinguishing between the two kinds of divergencies. The MS- 
scheme will be adopted for renormalization. 

The only difference between the non-Abelian and the Abelian theory on 
tree level is the color factor Cp = Tp{N 2 — 1)/N which multiplies the cou- 
pling constant in the potential, where N is the number of colors and Tp the 




(8) 




V(R) = - Urn — ln(Tr Pexp (ig jdx^ A^ a T a )) 



(9) 
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normalization of the generators, Tr[T a T b ] = 5 ab Tp. It is thus convenient to 
define 

= _ Cf ^M1, (10) 

q 

as this allows for an immediate generalization to sources in the adjoint rep- 
resentation: replacing Tfa — > —if akl and hence Cp —>■ Ca = N, the function 
ay describes the potential for static gluinos as well. 

On the one-loop level the difference between QED and QCD is more 
prominent. An obvious point is that the trilinear gluon self-coupling leads 
to a correction to the gluon propagator even if n/ = 0, and in principle to an 
additional vertex correction as well. But as a consequence of Feynman gauge 
and the special kinematics, every diagram containing a three gluon vertex 
with all three ends directly attached to the sources vanishes: if we denote the 
three-momenta flowing into the vertex with p, q, r, such a diagram involves 

v,v u v p (g^(p - qy + g^(q - r)" + g Pfi (r - p) v ) = 0. 

The same statement holds for the four-gluon vertex, which, however, first 
enters at the two-loop level. 

A second and more interesting point is that the color factors associated 
with the individual diagrams are not the same. Consider, for example, the 
ladder diagrams of Fig. again, 

Fig. |(a) ocCl , Fig. 0(6) oc C 2 F - C F ^. 

We can immediately identify the terms cx Cp as iterations of the tree-level 
potential, but there remains a term oc CpC& from the crossed ladder which, 
together with a corresponding term from the vertex correction that renders it 
infrared finite, leads to an additional contribution to the one-loop potential. 
This, of course, influences the way the exponentiation works at the two-loop 
level. 

As has been demonstrated in [[|, the consequence is that in order to 
compute the actual two-loop contribution to the potential, only those dia- 
grams have to be considered which involve color factors different from Cp 
and C f Ca and thus cannot result from iterations of the lower order diagrams. 
This means that for example the first three of the ladder diagrams in Fig. 
|^ and all graphs which are merely source self-energy insertions in one-loop 
graphs are irrelevant. 
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Figure 2: Two-loop ladder diagrams. 



To be more specific, the following diagrams have to be calculated: 

• The two-loop ladder diagrams Fig. |2|(d-f). 

• The corrections to the one-loop ladder diagrams shown in Fig. |3|. In 
general there would be more graphs of this type containing the three 
gluon vertex, e.g., analogous to Fig. 0(c), which, however, vanish in 
Feynman gauge as already explained. In fact, Fig. |3|(f) vanishes as 
well, but this is a consequence of considering the color-singlet state of 
the sources. 

• Two-loop vertex and gluon self-energy corrections, where the number of 
diagrams is also reduced by our choice of gauge, and double insertions 
of the corresponding one-loop corrections. 

• The graphs containing the four-gluon vertex with all ends attached to 
the sources vanish in Feynman gauge as well. 

All other two-loop graphs are already accounted for by the exponentiation. 

The relevant diagrams can be evaluated in momentum space without en- 
countering any special difficulties. Using the integration by parts method ||, 
most of the integrals that occur can be reduced to products or convolutions 
of the standard one-loop two-point function, its HQET equivalent as given 
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Figure 3: Corrections to one-loop ladder diagrams required for the potential 
calculation. 

in [|l(J, the HQET three-point function as given in [jllj and the mixed-type 
three-point function 

which can be computed by standard methods for vp = 0. The diagrams 
|2|(d-f), ^(e) and the vertex correction containing two three-gluon vertices, 
however, require the computation of some true two-loop integrals. As a 
detailed description of the calculation must be postponed to a future pub- 
lication, we only mention that the computer program FORM has been 
used for the evaluation of most of the diagrams and immediately present the 
results. 

Combining the two-loop result with the tree-level and one-loop expres- 
sions, the effective coupling constant introduced above can be written as 

^_„, , 2xl\ , «Ms(q 2 ) 31C A - 20T F n f 



MS 



4tt 9 



+( 



a^q 2 ) x 2 r _ 2 / 4343 . 24vr 2 - vr 4 . 22 



MS 



47T 



° A yi62 + i + y Cs 

1798 56 



-C A T F n f (— + yCa 



-C F T F n f (^-16( 3 ) + (T F n f f^ 



(11) 



By inserting the three-loop running coupling in the MS-scheme (the formula 
can be found, for example, in |TB|), we thus obtain the three-loop potential, 
in the sense that the expression is correct up to a constant multiplying a 3 . 
The terms proportional to C F T F and T F in (|ITD could have been obtained 
from the one- and two-loop gluon propagator, but the other two terms really 
required computing. 

Eq. (|TT|) can be used to determine the scheme-dependent coefficient fa of 
the /3-f unction for the V^-scheme, as defined by 



1 da 



v 



ay d In /i 2 



n=0 



(12) 



with the result (the first two coefficients, of course, coincide with those of 
the MS-scheme) 
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C\T F n f 



C A {T F n f f- 
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+2C 2 F T F n f + 



184 - 192C 3 
9^ 



C F (T F n f ) 2 . 



(13) 



The relation between the two couplings can, of course, be inverted easily, 
yielding 



a (q 

MS 



av(q 



a v (q 2 ) 31C A - 20T F n f 



«v(,q 

47T 



+C A T F n f 



r 2 



47T 

499 24tt 2 
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,400 



(14) 
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This formula could be used to improve the precision when extracting ql_(M§) 



from measurements of the Wilson loop on the lattice [pL4 
For N = 3 the results in actual numbers read 



a (q 2> 

T\,TQ V * 



a v (q 2 ) = «_(q 2 )|l + -^-^( 2 .583 - 0.278n/) 

+( ^ q (39.650 - 4.14771/ + 0.077nj)|, (15) 



7T 

2-, , 2\li «v(q 2> 



a_(q^) = «v(q z )<| 1 - -^-^(2.583 - 0.278n/ 



a v (q 2 )\2 



7T 



26.303 - 1.27772/ - 0.077n}j \, (16) 



which shows that the coefficients of the second order terms are not small 
even for rif = 5. 
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and for carefully reading the manuscript. This work was supported by the 
"Landesgraduiertenforderung" at the University of Karlsruhe. 
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